#include "../../include/ivpp.h"
#include <iostream>
#include <iomanip>
#include "analitic.h"
#include "pk_ivp.h"
#include <UTIL/TIMER/pentime.h>

using namespace std;

int main()
{
 float te=10;
 Tpoint_kin ivp;
 //Texplicit_euler<Tpoint_kin_state> ode(&ivp);
// Timplicit_euler<Tpoint_kin_state> ode(&ivp);
 Timplicit_tr<Tpoint_kin_state> ode(&ivp);
 //Trk4<Tpoint_kin_state> rk4ode(&ivp);
 //Tbdf<Tpoint_kin_state> bdfode(&ivp);
 *((Tpoint_kin_params*)&ivp)=TR;
 ode.dt=0.1;
 ode.t=0;
 ivp.x0(1.,ode.X());
 ivp.RO=0.1*ivp.betta;
  pentimer pt;
  pt.reset();
  pt.start();
 for(;ode.t<te;ode++)
 {
//  cout<<setw(5)<<ode.t<<' '<<ode.X()<<endl;
 }
  pt.stop();
  printf("%ld\n",pt.tim);
 return 0;
}

